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Abstract 

Arsenic, though a poor mutagen, is an accepted environmental carcinogen. Perturbation of DNA methylation 
pattern leading to aberrant gene expression has been hypothesized as the mechanism for arsenic induced 
carcinogenesis. We had earlier demonstrated the hypermethylation of promoter region of p53 and pi 6 genes in 
persons exposed to different doses of arsenic. Till now no genomic hot spot has been identified which is frequently 
hypermethylated or hypomethylated in persons chronically exposed to environmental arsenic. In the present work, 
we have identified one hypermethylated sequence by methyl-sensitive arbitrarily primed polymerase chain reaction 
in the peripheral blood leukocyte DNA of chronically arsenic exposed persons with and without arsenic induced 
skin cancer. The sequence is from GMDS gene responsible for fucose metabolism. Southern hybridization of the 
sequence to the amplification products of methyl sensitive restriction enzyme digested genome of persons 
exposed to different doses of arsenic indicated that methylation increased in a dose dependent manner. 
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Introduction 

According to (International agency for research on 
cancer 1997) and National Research Council (NRC 1999) 
arsenic is an important environmental toxicant and 
carcinogen. However, the mechanism of arsenic mediated 
carcinogenesis is not clear as arsenic is a poor mutagen 
(Rossman et al. 1980; Jacobson and Moltanbano 1985; 
Lee et al. 1985) and does not induce significant point 
mutations. Biotransformation of arsenic, on the other 
hand, involves methylation of inorganic arsenic to 
organic monomethyl arsonic acid (MMA) and dimethyl arsi- 
nic acid (DMA), using the same methyl donor S-Adenosyl 
methionine (SAM) also involved in DNA methylation 
(Vahter 1999). The interference of the DNA methylation 
pathway with arsenic detoxification pathway, as both 
the pathways require SAM, can lead to aberrant DNA 
methylation, resulting aberrant expression and/or silen- 
cing of genes Goering et al. (1999). Therefore, epigenetic 
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alterations, particularly aberrant DNA methylation has 
been mooted as a possible mechanism of arsenic induced 
carcinogenesis (Ren et al. 2010; Reichard and Puga 2010). 

Cytosine-5 methylation at the CpG islands in the regula- 
tory sequence of a gene is one of the key mechanisms of 
gene inactivation. DNA methylation/demethylation seems 
to regulate a plethora of biological processes involving 
transcription, differentiation, development, DNA repair, 
recombination, and chromosome organization. Perturbation 
of DNA methylation has been correlated with many cases 
of cancer (Jones and Baylin 2002). The hypothesis that 
arsenic perturbs DNA methylation has been tested success- 
fully on tissue culture system (Mass and Wang 1997), and 
later we demonstrated hypermethylation of the promoter 
region of p53 and pl6 genes in DNA extracted from 
peripheral blood leucocytes of persons exposed to 
different doses of arsenic (Chanda et al. 2006). A few 
highly exposed persons also showed p53 hypomethylation 
(Chanda et al. 2006). Further arsenic induced genome wide 
hypermethylation has been demonstrated by us in DNA 
extracted from same population Majumder et al. (2010). 

In this report we have further evaluated the hypothesis 
on a subsection of exposed population studied by isolating 
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hypermethylated sequence from genomic DNA of arsenic 
exposed persons by methyl sensitive arbitrarily primed 
polymerase chain reaction (MS-AP-PCR). Differentially 
methylalated fragments have been identified and isolated 
from chronically arsenic exposed people. 7 persons of ar- 
senic induced skin cancer out of 16 (all the cancer patients 
were recruited from previous study, Chanda et al. 2006), 
have one common hypermethylated fragment of 565 bp 
(this was cloned and sequenced). The same sequence was 
also isolated from 3 chronically arsenic exposed persons 
out of 10 (all of these 10 subjects were recruited from pre- 
vious study, Chanda et al. 2006). The sequence is then 
analysed by bioinformatic tools (NCBI BLAST) to indicate 
that the fragment is actually situated in the human GDP 
mannose 4-6 dehydratase gene (GMDS gene). Southern 
hybridization of this fragment to amplified products from 
methyl sensitive restriction enzyme digested genomic 
DNA of persons exposed to arsenic in drinking water indi- 
cated that the sequence is indeed hypermethylated. The 
product of the identified gene is involved in fucose metab- 
olism and it is reported that deletion of this gene results 
in cancer progression (Thompson et al. 1992; Becker and 
Lowe 2003; Yuan et al. 2008). Though have been found 
here in small proportion this hypermethylated fragment 
may be act as a potential target (probe) for detecting aber- 
rant methylation in chronic high level of arsenic exposure. 

Materials and methods 

Subject selection 

Subjects of this study were the same set of our earlier study 
on arsenic induced DNA hypermethylation in p53 and pl6 
gene promoter region and are all residents of South & North 
24 Parganas, West Bengal, India (Chanda et al. 2006). 

Criteria of diagnosis of arsenicosis and its severity are 
based on the parameters described earlier (GuhaMazumder 
et al. 1998; GuhaMazumder 2001; Chanda et al. 2006). In 
this study only the subjects for p53 gene promoter 
hypermethylation group of the previous study had 
been chosen. Participants had been divided into the 
five groups A, B, C, D according to the concentration 
of arsenic in their drinking water, i.e. 0-50, 51-250, 
251-500, 501-1000 ug/1 respectively as earlier and 
group E with 500-1100 ug/1 of arsenic suffering from 
arsenic induced skin cancer. As the concentration of 
arsenic in group A is within the permissible limit according 
to WHO and Medical Council of India, it was considered 
as the unexposed control group (NRC 1999). Initially the 
number of participants in each group was 24, 12, 18, 
15 and 16 in A, B, C, D and E group respectively 
(Chanda et al. 2006). Among those, 11 subjects in group C, 
10 subjects in group D and all the 16 subjects in group E 
were chosen for this study. All of the subjects 
chosen for MS-AP PCR have hypermethylated p53 
promoter region compared to normal unexposed 



persons (Chanda et al. 2006). All of the subjects studied for 
MS-AP-PCR were compared to normal unexposed subjects 
treated in similar way. In this study 12 of the normal 
unexposed subjects were recruited from group A. 

The initial studies on isolation of hyper/hypo methylated 
stretch of genomic DNA was performed with peripheral 
blood leukocyte DNA of highly arsenic exposed persons of 
group D and arsenic induced cancer group, E. Later, lower 
exposure group C was also evaluated for the presence of 
such hyper/hypomethylated gene fragments. Among 16 of 
the cancer patients (group E) studied, 7 had a hypermethy- 
lated DNA fragment of 565 nucleotide long sequence. 
Among 10 subjects of group D, 3 had that hypermethylated 
fragment. The fragment identified was from the partici- 
pants of group D and group E but not from any lower 
exposure group. Although there is an overlap between 
group D and E in respect to the concentration of arsenic in 
water but the difference is one group have arsenic induced 
skin cancer with higher degree of skin manifestations 
(group E) while the other group (group D) is only character- 
ized by higher degree of skin manifestations without cancer. 
The identified and isolated hypermethylated fragment was 
then sequenced. The sequence is from the intronic region 
of human GMDS gene situated in between exon 1 and 2. 
The southern hybridization studies of the identified frag- 
ment with DNA of 4 persons, taken 1 from each exposure 
group indicate that the sequence is indeed hypermethylated. 
Demographic data for this study population is described in 
detail (Table 1). Once the fragment was identified and 
isolated from peripheral blood leukocyte DNA of arsenic 
induced cancer patients, the procedure was cross-checked 
using DNA samples isolated from cancer biopsy samples 
of the same patients. But it was not done in case of 
group D samples due to lack of biopsy tissues in those 
cases (as these are not arsenic induced cancer). 

Written informed consent was obtained from all 
participants before drawing their blood. The name of 
the institute where human clinical studies were carried 
out is Institute of Post Graduate Medical Education and 
Research, Kolkata (IPGME&R), which is run by Govt, of 
West Bengal, a state government within the framework of 
Republic of India. 

Molecular Biological and in silico experiments were 
carried out in University of Calcutta and Presidency 
University, Kolkata which are also run by Govt, of West 
Bengal. Ethical principles followed by the institute are 
guided by rules as formulated by Indian Council of 
Medical Research and these are in agreement with 
Helsinki declaration. 

Determination of Arsenic concentration in urine 
and water 

Level of arsenic in drinking water and urine was de- 
termined by atomic absorption spectrophotometer 
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Table 1 Demographic data of study subjects taken from different arsenic exposure groups 



Age group 


Sex 


Group A (0-50 ug/l) 
p53 methylation 


Group C (250-500 ug/l) 
p53 methylation 


Group D (501-1000 ug/l) 
p53 methylation 


Group E (500-1100 ug/l) 
p53 methylation 


<20 years 


Male 


N = 2; 0.08, 0.19 


N = 


= 2; 1 .34, 1 .62 








Female 




N = 


= 2; 1 .50, 1 .45 






21-40 years 


Male 


N = 3; 0.03, 0.23, 0.15 


N = 


= 3; 1.50, 2.78, 5.00 


N = 4; 1 .95, 2.04, 4.3 2.46 


N = 3; 0.85, 1.09, 1.32 




Female 


N = l; 0.15 






N = 1;2.56 


N = 3; 1 .00, 1 .62, 1 .09 


41-60 years 


Male 


N = 3,0.1 8, 0.27, 0.27 


N = 


= 2; 1 .56, 2.4 


N = 2; 2.85, 4.46 


N = 4; 1 .42, 3.08, 2.09, 2.09 




Female 


N = 3; 0.08, 0.09, 0.32 


N = 


= 2; 1 .95, 4.66 


N = 2 4.45, 2.7 


N = 1; 2.55 


>60 years 


Male 








N = 1; 2.15 


N=4; 2.20, 2.23, 1.62, 1.60 




Female 










N = 1;2.11 


Smoking status 


Smoker 


8 


9 




7 


11 




Nonsmoker 


4 


2 




3 


5 




Exsmoker 












Avarage duration of exposure 




1 1 .5 years 


15 


years 


10 years 


1 7 years 


Total number of samples 


47/M 


12 


11 




10 


16 



Note: numerical values in each cell indicate the degree of p53 methylation for individual study subjects recruited in the study. 



with hydride generation system (AAS) Atallah and 
Kalman (1991). 

DNA isolation from blood 

DNA was extracted from whole blood by conventional 
chloroform extraction method using 0.01% SDS and 
Proteinase K (0.1 mg/ml) (Miller et al. 1998). 

DNA isolation from tissue 

DNA was extracted from cancer biopsy tissue samples 
by conventional phenol-chloroform (1: 1, v/v) extraction 
method and then by chloroform extraction followed by 
salting out using 0.01% SDS and Proteinase K (0.1 mg/ml) 
(Miller et al. 1998). 

p53 methylation status analysis 

The p53 tumor suppressor gene methylation status was 
analyzed in each subject by the method described earlier 
(Chanda et al. 2006). 

Determination of clinical symptom score 

Each subject was assigned a clinical symptom score 
which reflects the severity of his/her skin manifestations. 
Both pigmentation and keratosis were graded as 1,2 or 3, 
depending on the level and severity of symptoms. Sum of 
the two was clinical symptom score, so that a person can 
have maximum score of 6. Control subjects have no 
pigmentation and keratosis and therefore have a clinical 
symptom score 0. The detail structure of the scoring 
system for pigmentation and keratosis is given in Table 2. 

Restriction enzyme digestion for arbitrarily primed PCR 

Concentration and quality of isolated genomic DNA was 
determined UV-vis spectrophotometer (OD 260/280 >1.8). 



300 ng of total genomic DNA isolated from persons, 
unexposed/exposed to arsenic through drinking water, 
was digested with 5 units of Rsal and 5units of Hpall 
restriction enzyme at 37°C overnight. Hpall is a methyla- 
tion sensitive isoshizomer of Mspl whose recognition 
sequence is CCGG. A sequence hypermethylated at this site 
would not be digested, whereas the unmethylated DNA 
would. The persons taken for MS-AP-PCR were from 
higher exposure groups of arsenic (251-500 ug/l, i.e. group 
C; 500-1000 ug/l, i.e. group D; and arsenic induced cancer 
group, E, with an exposure level of 500-1100 ug/l) and all 
have hypermethylated p53 promoter. Out of 18 in group C, 
11 subjects were taken for MS-AP-PCR. All have hyper- 
methylated p53 promoter region. In group D only 10 sam- 
ples were chosen with hypermethylated p53 promoter 
having a median value of 2.63. In group E all the 16 samples 
were studied with pS3 promoter hypermethylation 
with a median value of 1.62. The median value for p53 
methylation in group A (unexposed control group) was 
0.26 which is treated here as a basal value for normal 
unexposed persons (Chanda et al. 2006). Demographic 
data and p53 methylation values for subjects included in 
this study are presented in Table 1. 

Methyl sensitive arbitrarily primed PCR (MS-AP- PCR) 

When Rsal + Hpall digested DNA was used as template 
in MS-AP-PCR using random primers that target 
CG-rich DNA sequences (Zhong and Mass 2001), a 
series of amplified products were observed. Of these, 
a band present in PCR products of arsenic exposed DNA 
but absent in PCR products of similarly digested unexposed 
DNA represents the region of hypermethylation (Zhong 
and Mass 2001). We used 3 different primers. Amongst 
these primers, OPN Hindl2 (5'-AGCTTCTCCCTC-3') 
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Table 2 Dermatological criteria and graduation of chronic arsenic toxicity for scoring sustem of skin manifestations 

Pigmentation 



Mild 1 

Defuse Melanosis, Mild Spotty pigmentation, 
Leucomelanosis 

Keratosis 

Mild Score = 1 

Slight thickening, or minute papules (<2 cm) 
in palm and soles 



Moderate score = 2 

Moderate Spotty pigmentation 

Moderate score = 2 

Multiple raised keratosis papules (2 to 5 cm) 
in palm & soles with diffuse thickening 



Severe score = 3 

Blotchy Pigmentation, Pigmentation of under 
surface of tongue, buccal mucosa 



Severe score = 3 

Diffuse severe thickening, large discreet or 
confluent keratotic elevations (>5 cm), palm 
and soles (also dorsum of extremely and trunk) 



The underlined data represents the clinical symptom score. 



gave one common hypermethylated fragment in arsenic 
induced cancer patients and in highly arsenic exposed 
persons when subjected to PCR amplification. The 
concentration of OPN Hind 12 in the PCR reaction 
mixture was 0.5uM. The PCR protocol was initial de- 
naturation at 94°C for 5 minutes, 35 cycles at 94°C 
for 1 min, 40°C for 1 min, 72°C for 2 min, followed 
by 10 min at 72°C. 

Isolation of candidate bands 

The PCR products from DNA of arsenic exposed persons 
were compared with PCR products from control DNA. 
The band, which appears only in the exposed DNA but 
not in unexposed DNA was supposed to be region of 
hypermethylation. Similarly, band which appears only in 
the unexposed control but not in the exposed DNA indi- 
cated the site of hypomethylation in the DNA of exposed 
persons. Such a region of hypermethylation from chronic- 
ally high arsenic exposed people with and without arsenic 
induced cancer was identified. The ethidium bromide 
stained PCR amplified DNA band was then excised from 
the gel by a scalpel and recovered by the usual 'crush and 
soak' method (Sambrook et al. 1989). The candidate band 
isolated from subjects were from group D and E. Clinical 
symptom score, pS3 methylation status and degree of 
arsenic exposure for those subjects are given in Table 3. 

DNA cloning in plasmid vector 

The gel- recovered PCR product was re-amplified using 
the same PCR protocol with same primer. The amplified 
product was then purified by ethanol precipitation and 
cloned in E.coli XL1 blue strain using pTZ57R/T vector 
(TA cloning kit, Fermentas). The positive clones were 
identified by performing colony PCR with universal 
primers and sequenced. 

Southern hybridization 

Exacdy equal amount (1.3 ug) of genomic DNA of four per- 
sons from four different exposure groups (Group A, B, C, D) 
were subjected to restriction digestion by Rsal and 
Hpall and incubated overnight at 37°C. Each of the 



digested products was then subjected to PCR amplification 
using primer OPN Hindl2. The PCR products obtained 
from four different DNA samples were resolved by electro- 
phoresis on a 2% agarose gel. The gel was blotted on nylon 
membrane using standard technique and then hybridized 
with a- P 32 dCTP (BARC, India) labelled clone insert. The 
same procedure of hybridization was carried out using one 
DNA sample from cancer patient where instead of group 
A, B, C, D group B,C, D and E were used to hybridise with 
the labelled clone of the fragment isolated. The relevant 
parameters for the persons taken from four different 
groups for hybridization are described in Table 4. 

Result 

Using the technique of MS-AP PCR, 1 common hyper- 
methylated DNA fragment was identified from 10 different 
people with chronic high level of arsenic exposure with and 
with out cancer. Among 16 of the arsenic induced cancer 
patients studied (belonging to group E) 7 have the hyper- 
methylated DNA fragment of 565 nucleotide base pair. 
Among 10 of group D subjects 3 have been identified to 
harbour this hypermethylated DNA fragment. Demographic 
data and p53 methylation status for these 10 subjects 
(with hypermethylated DNA) has been listed in Table 3. 
Interestingly, people from lower arsenic exposure (group C) 
did not have this hypermethylated gene fragment. 

The fragment identified is a region of hypermethylation 
in comparison to normal unexposed persons. DNA 
sequence analysis revealed that the identified fragment 
has significant homology match (99%) to the sequence 
of human GMDS gene (Accession no. NT_007592.15), 
Homo sapiens) (taken from GENEBANK database) 
after BLAST search. The sequence is situated in the intron 
between exon 1 and 2 of GMDS gene. (Genomic context: 
chromosome: 6; Maps: 6p24.1-25.3). It is the longest 
intronic sequence in GMDS gene (> 1,80,000 bp). This gene 
is involved in carbohydrate metabolism and generation of 
fucose. Fucose mediates initial contact between extravagat- 
ing leucocytes and endothelial cells. Influence of fucose gen- 
erating enzymes on leukocyte adhesion activity has been 
reported (Sullivan et al. 1998; Eshel et al. 2001). 
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Table 3 Demographic data and p53 methylation status of subjects having GMDS gene hypermethylation 

Sample Age (yr)/sex Smoking status Cone, of arsenic Duration of Degree of Total urinary p53 methylation value 

in water ug/l exposure yrs pigmentation arsenic ug/l 



KA 261 


52/M 


Non smoker 


580 


7 


++ + 3 


272.8 


4.46 


Gr. D 
















DHW 088 


43/M 


smoker 


683 


5 


++++ 4 


89 


2.85 


Gr. D 
















CW045 


38/M 


smoker 


531 


5 


++++ 4 


189 


2.46 


Gr.D 
















CNBB 33 


48/F 


Non-smoker 


826 




++++ 4 


212 


2.55 


Gr.E 
















CNBB 28 


40/M 


Ex smoker 


740 


10 


++++ 4 


126 


1.62 


Gr.E 
















A 10 


51/M 


smoker 


514 


1/ 


++++ ++ 6 


211 


3.08 


Gr.E 
















A 15 


47/M 


smoker 


623 


13 


++++ ++ 6 


97 


2.09 


Gr.E 
















A 20 


53/M 


smoker 


744 


10 


++++ 4 


143 


2.09 


Gr.E 
















A 17 


63/M 


smoker 


556 


17 


++++ ++ 6 


171 


2.20 


Gr.E 
















A 21 


61 M 


smoker 


631 


10 


++++ ++ 6 


206 


2.23 


Gr.E 

















Note: The pigmentation and keratosis was assigned as a numerical score according to the degree of severity. 



When the insert of the clone is hybridized to the PCR 
products amplified by the OPN Hind 12 primer from Hpall 
digested genomic DNA of persons exposed to various doses 
of arsenic and arsenic induced cancer, it was found that the 
hybridization increases in higher exposure groups and in 
arsenic induced cancer patients. This indicates that the 
segment is indeed hypermethylated in genomic DNA of 
persons with high arsenic exposure and also in arsenic 
induced cancer patients (Figure la, lb). Hypermethylation 
rendered the fragment insensitive to digestion by the methyl 
sensitive enzyme Hpall at the relevant site in higher expos- 
ure group DNA and the desired region was available for 
amplification. So the amount of PCR product template avail- 
able for associating with the probe is more in the arsenic 
exposed group and in cancer group with hypermethylation 
in their p53 promoter region than in the unexposed group. 



Our identified fragment, OPN Aga 8, shows very low 
association with the DNA of <50 ug/l exposure group 
(Figure la). The degree of association increases gradually 
with the degree of arsenic exposure in higher exposure 
groups. The fragment isolated has a higher degree of 
association during hybridization with PCR product of 
person with arsenic induced cancer (Figure lb). 

Our results indicated that in the GMDS mRNA dataset, 
non-linearity was found when pair wise sequence diver- 
gence per site corrected from Kimura-2-parameter (K2P) 
distances were plotted against p-distances for both transi- 
tion and transversion among the studied animal species, 
but Ts sites exhibited slightly higher saturation over 
Tv (Figure 2a and 2b). We have estimated Disparity 
Index (DI) per site of GMDS mRNA sequences between 
all sequence pairs Kumar and Gadagkar (2001). Value 



Table 4 Demographic data of subjects taken from different exposure group for southern hybridization 


Sample Age (yr)/ 
sex 


Smoking 
status 


Cone, of arsenic in 
water (ug/l) 


Duration of 
exposure yrs 


Degree of 
pigmentation 


Total urinary 
arsenic (ug/l) 


p53 methylation 
value 


1 40/M 


smoker 


11 


5 


-(0) 


0.0 


0.20 


2 47/M 


smoker 


118 


5 


+ (D 


31 


0.78 


3 52/M 


smoker 


314 


7 


+++ (4) 


87 


1.95 


4 39/M 


smoker 


644 


6 


++++ (6) 


223 


2.56 
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>500 251- 51- <50 
500 250 




OPN Aaa8 



B 

50 -250 



C 

250 - 500 



D 

500 -1000 



E 

500 - 1100 




Figure 1 Represents the southern blots of cloned insert, a. Southern hybridization pattern of the cloned insert (OPN Aga8) to amplification 
products of Hpall digested DNA from four persons of different exposure groups, b. Southern hybridization pattern of the cloned insert (OPN Aga 8) 
to amplification products of Hpall digested DNA from four persons of three different exposure groups without arsenic induced cancer and one group 
of arsenic induced cancer. 



greater than '0' indicate larger difference in base compos- 
ition biases, than expected. This is based on evolutionary 
divergence between sequences and by chance alone. There 
were a total of 670 positions in the final dataset. Highest 
Disparity Index was observed between Apis florae 
and Anopheles gambiae (DI = 24.38) among the other 
sequence pairs (Figure 3). 

Phylogenetic analysis 

Human GMDS gene sequence comparison was done with 
42 different species of different animal groups, for which 
the sequences available in GenBank gives us the evolution- 
ary relationship of this gene among the selected species 
(Table 5). Multiple sequence alignment was executed with 
three dataset derived from GMDS mRNA sequences using 
Clustal W and it was found that the sequence identified is 
conserved in a number of genera studied. The sequence 



was further analyzed using the Kimura 2-parameter model, 
p-distance to assay the probability of the number of transi- 
tional and transversional substitutions per site between 
sequences. All positions containing gaps and missing 
data were eliminated. Phylogenetic tree was constructed 
based on Neighbor-Joining method (NJ) with Kimura 
2-parameter using MEGA version 5.05 Tamura et al. 
(2011); Saha et al. 2013a, b) from both transition and trans- 
version data. Standard error estimate(s) were obtained by 
bootstrap procedure (1000 replicates). Disparity Index 
per site was estimated for all sequence pairs Kumar and 
Gadagkar (2001). 

Maximum Composite Likelihood Estimate of the 
pattern of nucleotide substitution was estimated according 
to Tamura et al. (2004) where each entry shows the prob- 
ability of substitution (r) from one base (row) to another 
base (column) (Table 6). For simplicity, the sum of r values 
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0.000 0.100 0.200 0.300 0.400 0.500 0.000 0.050 0.100 0.150 0.200 0.250 0.300 

K2P Distance K2P Distance 

Figure 2 Pairwise sequence divergence among the 42 animal taxa. GMDS mDNA, (a) and (b); plots of Kimura 2 parameter (K2P) inferred 
Transition (Ts) and Transversion (TV) distances against the P-distance. 



is made equal to 100. Rates of different transitional 
substitutions are shown in bold and those of trans- 
versional substitutions are shown in italics. The nucleo- 
tide frequencies are 28.32% (A), 26.46% (T/U), 24.03% (C), 
and 21.20% (G). The transition/transversion rate ratios are 
kl = 1.887 (purines) and k2 = 3.132 (pyrimidines). The 
overall transition/transversion bias is R =1.219, where 
R = [A*G*« + T*C*£2]/[(A + G)*(T + C)]. The whole analysis 
involved 42 different nucleotide sequences. All positions 
containing gaps and missing data were eliminated. There 
were a total of 670 positions in the final dataset. 

GMDS mRNA derived phylogenetic tree inferred by 
the NJ method represented in Figure 4. Due to saturation 
of both the substitutions, the sum of the transition and 
transversion for phylogenetic tree reconstruction by the 
NJ method based on K2P model has been used. These 
sequences are composed of 435 variable sites and 389 
parsimony informative sites. The transition/transversion 
ratio and the overall mean distance has been found to be 



0.96 and 0.333 ±0.016. We have observed that GMDS 
sequence of Callithrix jacchus shared a common 
ancestor with the sister clade containing Homo sapiens, 
Pan troglodytes, Pongo abelii, Nomascus leucogenys and 
Macaca mulatta with 100% bootstrap support. Equus 
caballus have shown monophyly with closely related sister 
species Ailuropoda melanoleuca and Canis lupus which 
was supported by high bootstrap value (97%). Cricetulus 
griseus, Rattus norvegicus and Mus musculus formed a 
monophyletic group with high bootstrap support 
(100% and 99% respectively), whereas Ornithorhynchus 
anatinus diverged early in the tree among all other 
mammamls under study. NJ tree also depicted that 
Taeniopygia guttata, Gallus gallus and Meleagris gallopavo 
belong to the class Aves that have exhibited mono- 
phyletic origin and evolved parallely with the reptiles 
(Anolis carolinensis) but diverged after the class Amphibia 
and Actinopterygii. Class Insecta belongs to the phylum 
Arthropoda consisted of two clades, one of order Diptera 




lOoiONT-inoinsT-ioronsi-inoinNi-iooinsi- 
co(OOcor^ocoi^o-*r^O'*i^i-'^-i^'i-'*oo-i-'a-coi-m 
i-T-i-Nwwncon^**in«)inoiB(DSSMOco 



Sequence pair 

Figure 3 Disparity Index per site is shown for all sequence pairs for GMDS sequences. All positions containing gaps and missing data were 
eliminated. Values greater than 0 indicate the larger differences in base composition biases than expected based on evolutionary divergence 
between sequences and by chance alone. The analysis involved 42 nucleotide sequences. There were a total of 670 positions in the final dataset. 
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Table 5 GenBank accession numbers and size of GMD 



sequence of sampled taxa 



Sr. No. 


Organism 


Accession No. 


Size 


1 


Homo sapiens 


BC000 117.1 


1119 bp 


2 


Pan troglodytes 


XM_5 18203.3 


795 bp 


3 


Pongo abelii 


XM_0028 16345.1 


1119 bp 


4 


Nomascus leucogenys 


XM_003272184.1 


1119 bp 


5 


Macaca mulatto 


NM_00 1266789.1 


1119 bp 


6 


Callithrix jacchus 


XM_002746279.2 


1119 bp 


/ 


Equus caballus 


XM_001 490703.3 


1050 bp 


8 


Ailuropoda melanoleuca 


XM_002922834.1 


1119 bp 


9 


Canis lupus 


XM_545311.3 


1011 bp 


10 


Loxodonta africana 


XM_00341 7823.1 


1119 bp 


11 


Cavia porcellus 


XM_003463230.1 


1050 bp 


12 


Bos taurus 


NM_001 080331.1 


1119 bp 


13 


Oryctolagus cuniculus 


XM_002720970.1 


1662 bp 


14 


Cricetulus griseus 


NM_00 1246696.1 


1119 bp 


15 


Rattus norvegicus 


NM_00 1039606.1 


1119 bp 


16 


Mus musculus 


BC093502.1 


1119 bp 


17 


Ornithorhynchus anatinus 


XM_001 51 0089.1 


1260 bp 


18 


Anolis carolinensis 


XM_003225586.1 


1006 bp 


19 


Taenlopygia guttata 


XM_0021 97547.1 


1 035 bp 


20 


Gallus gallus 


XM_41 8977.3 


1086 bp 


21 


Meleaghs gallopavo 


XM_003204780.1 


1047 bp 


22 


Xenopus laevis 


BC1 5741 1.1 


1110 bp 


23 


Danio rerio 


NM_001 102475.2 


1113 bp 


24 


Salmo salar 


NM_001 141373.1 


1113 bp 


25 


Oreochromis niloticus 


XM_00345 7295.1 


1116 bp 


26 


Nematostella vectensis 


XM_00 1622499.1 


1077 bp 


27 


Trlchoplax adhaerens 


XM_0021 161 09.1 


1080 bp 


28 


Brachionus manjavacas 


FJ829249.1 


1027 bp 


29 


Culex quinquefasciatus 


XM_00 1868832.1 


1107 bp 


30 


Anopheles gambiae 


XM_308963.3 


1089 bp 


31 


Aedes aegypti 


XM_001 650058.1 


1149 bp 


32 


Tribolium castaneum 


XM_968229.1 


1071 bp 


33 


Drosophila willlstoni 


XM_002066598.1 


1194 bp 


34 


Amphimedon queensiandica 


XM_003384374.1 


1110 bp 


35 


Brugia malayi 


XM_00 1898680.1 


1 1 64 bp 


36 


Loa loa 


XM_003138092.1 


1143 bp 


37 


Dictyostelium purpureum 


XM_003283436.1 


1068 bp 


38 


Acyrthosiphon pisum 


XM_00 1949034.2 


1 086 bp 


39 


Nasonia vitripennis 


XM_001 605356.2 


1071 bp 


40 


Megachlle rotundata 


XM_003 702009.1 


1077 bp 


41 


Bombus impatiens 


XM_003484683.1 


1071 bp 


42 


Apis florea 


XM_003692995.1 


1077 bp 



and the other of Hymenoptera. Two Nematod species, 
Brugia malayi and Loa loa are closely related providing 
90% sequence similarity. 

Identified sequence OPN Aga 8 

TCCCTCACTA CTCAAAGTTG ATGACTTCTT 

AAACCAAAAT GGTTGGTCAG AATCCAATCA 
AGAATATAAA GGCAACTGAA TAAATAAAAC CAT- 

AAAGTAA GTGTAAAATA CTAGTGTCCA 

GACATCTGAG ATGTATGTGG CTACTATGAA 

ACTTCCACAG CTGTACCGGC CGGGAGCTCA 

CGTGGTTCCC CAGGTTTAAC AGAACCCATT 

ACCAGTAAGA GTTTTATTTG CTTAATAAAT 

TACATTCTAA AG C AC AATAG CCTAGGCTCA 

TAGCTGTAAA ATTGCCAAAT ATTGTCAATG 

ACCACTCTCT GGTCATAAAT AACAAAATAA 

TCTTGTGACT CATTGGATTT TTGATTCC- 

CAAGGCGATTCT TTCTCGCCAT TACTCAAAAA 
TGTGAAAAAG TGCCTCTACGTGGCATTTTA 

TGGAGGATAT AAATTACTCA AAGGAGATGA 

CATAGGACAGATTTGTAGGC CGAGTAACAG 

GAACCAGCCA ACCAACTGTG TAAATTAAA- 

GAACTAGTGAC AAAGAAGAGG GCTAGTGAAA 
GAATTCTGAA ATCCTAAGAA CAGAT 

Discussion 

The mechanism by which arsenic contributes to the 
development of cancer is currently a subject of intense 
interest. Arsenic does not act as a point mutagen. 
However, metabolism of arsenic involves methylation 
of inorganic arsenate to dimethyl arsinic acid via alternating 
reduction of pentavalent arsenic to trivalent arsenic and 
addition of methyl group (Vahter 1999; Donohue and 
Abernathy 2001). The arsenic methyl transferase uses 
the same methyl donor SAM as DNA methyltransferase 
(Dnmt) and other methyltransferases. Interaction of arsenic 
methylation/detoxification pathway with DNA methylation 
pathway and consequent imbalance in DNA methylation 
has been envisaged. Increase of cytosine methyltransferase 
transcript after arsenic exposure has been reported 
(Zhong et al. 2001), and this might explain the initial 
hypermethylation through excessive induction of the 
enzyme. On the other hand prolonged arsenic exposure 
may cause depletion of the SAM pool due to over con- 
sumption of the methyl groups by arsenic methyltransfer- 
ase, and cause hypomethylation of DNA. Although we have 
failed to isolate any hypomethylated fragment from patients 
who have arsenic induced cancer or persons having chronic 
high level of exposure with systemic manifestations, yet, 
there was number of subjects with p53 promoter hypome- 
thylation in our previous study (Chanda et al. 2006). In fact 
decrease of tissue arsenic burden has been correlated with 
methionine intake in experiments with laboratory rats 
exposed to arsenic (Nandi et al. 2005). Interestingly the 
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Table 6 Maximum composite likelihood estimate of the 
pattern of nucleotide substitution 

A T C G 

A - 5.91 4.73 10.12 

T 6.32 - 14.82 5.36 

C 6.32 18.5 - 5.36 

G 11.93 5.91 4.73 

Note: Specificity for the bold symbols are justified in result. 

fragment of GMDS gene isolated from subjects having high 
degree of arsenic exposure and relatively high degree of p53 
methylation are from male subjects except one. We had 
not found any correlation between sex and pS3 methylation 
status in our previous study (Chanda et al. 2006). The 
hypermethylated gene fragment has been isolated from 
both smokers and nonsmokers although in the present 
study the number of smoker having GMDS gene intron 



hypermethylation is 8 out of 10. In our previous study we 
showed that there was no association between the p53 and 
pl6 methylation status with individual's smoking habit. In 
this study, due to small sample number, we could not carry 
out statistical evaluation for correlation between smoking 
habit and GMDS gene methylation. The fragment isolated 
in this study is from 9 male subjects and only one female 
subject out of 10 subjects. Increase in number of study 
subjects and consequent rise in the number of subjects 
having GMDS gene hypermethylation may provide 
definite information about the association (if any) of 
sex specificity of GMDS gene hypermethylation with 
arsenic exposure in human. 

Such aberrant methylation of the genome leading to 
gene expression anomalies and have been mooted as 
possible mechanism of arsenic induced carcinogenesis. 
Cancer, which results from inappropriate expression of 
genes, may involve hypomethylation of protooncogene 
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Figure 4 Unrooted neighbor-joining tree constructed from the GMDS mDNA sequences (branch length = 4.18962548). The percentage 

of replicate trees in which the associated taxa clustered together in the bootstrap test (1000 replicates) is shown above the branches. The 

evolutionary distances were computed using the Kimura 2-parameter (K2P) method and are in the units of the number of base differences per 

site. There were a total of 670 positions in the final dataset. 
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and/or hypermethylation of tumor suppressor genes that 
can alter the level of their expression and thereby promote 
cancer. In fact, there are recent observations that 
widespread methylation changes occur during tumor 
development (Jones and Baylin 2002). 

Overall, tumor cell DNA is hypomethylated compared 
to normal cell DNA and underexpression of Dnmtl gene 
causes aggressive tumor induction in genetically engineered 
mice (Gaudet et al. 2003). However, for some tumor 
suppressors like pl6, pi 5 methylation is a common 
alternative to point mutation and in others like RASSF1A 
or H1C1, it is the only mechanism for tumor specific loss 
of function (Jones and Baylin 2002). Silencing of genes like 
TIMP-3 through methylation has been associated with 
metastasis (Darnton et al. 2005). 

Methylation of DNA is maintained by a balance of the 
activity of DNA methyltransf erase (Dnmt 1, Dnmt 3a 
and Dnmt 3b) and DNA demethylase (mbd2) activity. 
Inhibition of mbd2 by antisense expression results in 
inhibition of anchorage-independent growth of antisense 
transfected cancer cells or cells infected with an adeno- 
viral vector expressing antisense mbd2 (Slack et al. 2002). 
Expression of Dnmt mRNA is significantly high in gastric 
cancer in comparison to non-cancerous gastric mucosa 
(Fang et al. 2004). Similarly level of mbd2 mRNA level is 
significantly lower in gastric cancer tissue than normal 
gastric mucosa (Fang et al. 2004). 

Previous works with human adenocarcinoma cell 
line in tissue culture showed that arsenic induces 
significant changes in methylation status in tumor 
suppressor gene p53 Mass and Wang (1997). Later, 
using arsenic exposed human kidney cell lines global 
hyper and hypomethylation has been demonstrated by 
the same group (Zhong et al. 2001). DNA sequencing 
and SssI methylase assay were used for estimation of 
genomic CpG methylation level. Arsenic exposure of 
A 549 cells in culture resulted in a dose dependent 
increase in cytosine methylation in p53 gene and a 
small increase in global methylation Mass and Wang 
(1997). Later we have shown that arsenic induces 
genomic hypermethylation in chronically exposed persons 
Majumder et al. (2010). An increase in the rate of tran- 
scription of DNA methyltransferase gene in cells exposed 
to arsenite was detected by RT-PCR (Zhong et al. 2001). 
Our group has demonstrated for the first time that there 
is dose dependent enhancement of methylation in the 
promoter region of p53 and pl6 tumor suppressor 
genes of genomic DNA extracted from peripheral 
blood leucocytes of persons exposed to various doses 
of arsenic (Chanda et al. 2006). However, both these genes 
are associated with cell cycling and repair, and the possibility 
exists that methylation perturbation observed is engineered 
through disturbances in cell cycle produced by arsenic, 
and is local, rather than a global effect of arsenic. 



In the present work we have investigated that whether 
there is any probable common target for aberrant DNA 
methylation after arsenic exposure in exposed persons 
apart from p53 or pl6 gene methylation. We have 
successfully identified one fragment of hypermethylated 
DNA from persons exposed chronically to arsenic and per- 
sons having arsenic cancer. The subjects have been chosen 
from our previous study population (Chanda et al. 2006) 
having hypermethylated p53 promoter region with chronic 
high level of arsenic exposure with and without arsenic 
induced cancer. Therefore persons having GMDS gene 
intron hypermethylation also have p53 promoter hyperme- 
thylation. Thus this study reflects an association between 
the p53 promoter hypermethylation with GMDS gene 
intron hypermethylation in chronic high level of arsenic 
exposed people. The fragment was isolated from both per- 
ipheral blood leukocyte DNA and from cancer biopsy tissue 
of persons having arsenic induced cancer. The hypermethy- 
lated DNA fragment is from GMDS gene responsible for 
fucose metabolism. GMDS is the binding partner of tankyr- 
ase which is needed to be associated for the first step of 
fucose biosysnthesis (Bisht et al. 2012). Oligosaccharides 
are involved in various aspects of life process including 
birth, differentiation, growth, inflammation, carcinogenesis, 
and cancer metastasis. Fucosylation is one of the most 
important oligosaccharide modifications in cancer. This 
type of glycomodification can be treated as a biomarker in 
cancer (Moriwaki et al. 2009; Miyoshi et al. 2012). 

Fucosylated alpha-fetoprotein (AFP) is widely used in 
the diagnosis of hepatocellular haptoglobin have also 
been found in sera of patients with various carcinomas 
(Miyoshi et al. 2012). Deletion mutation of the GMDS 
gene plays a pivotal role in fucosylation in human colon 
cancer. Loss of function mutation of this gene may lead 
to a virtually complete deficiency of cellular fucosylation, 
tumor progression and metastasis (Nakayama et al. 2013) 
and transfection of the wild-type GMDS into HCT116 
cells restored the cellular fucosylation. This type of 
GMDS mutation resulted in resistance to 77?A/L-induced 
apoptosis followed by escape from immune surveillance 
(Moriwaki et al. 2009; Haltiwanger 2009; Moriwaki et al. 
2011) and thus promote carcinogenesis. Further, epigen- 
etic regulation of fucosylation and TRAIL induced apop- 
tosis in conjunction to cancer had been studied by same 
group (Moriwaki et al. 2010). Although in the present 
study we have not shown any association with the level of 
fucose in patients with GMDS gene hypermethylation, but 
still GMDS gene fragment hypermethylation is associated 
with p53 hypermethylation with development of arsenic 
induced cancer (in group E) or severe skin manifestations 
(in group D) as a result of chronic high level of 
exposure. In the present study we have not work out the 
degree of correlation (if any,) between the GMDS intron 
hypermethylation and p53 promoter hypermethylation, 



Chanda et al. SpringerPlus 2013, 2:557 
http://www.springerplus.eom/content/2/1/557 



Page 11 of 12 



(as quantitative analysis of GMDS gene hypermethylation 
has not been studied here), although this study signifies 
the association between p53 promoter hypermethylation 
and GMDS gene hypermethylation after chronic arsenic 
exposure. 

The intronic fragment isolated in the present study 
showed hypermethylation in comparison to normal unex- 
posed subjects. This epigenetic modification may be involve 
in transcriptional modification and may modify the ultim- 
ate cellular level of GMDS enzyme. As it is reported that 
aberrant methylation of introns or intergenic regions can 
regulate non coding RNA function to modify the degree of 
transcription of a gene and the exonal expression is 
dependent over the local methylation status rather than the 
promoter region (Cheung et al. 2011). Consequences of 
intron methylation have also been studied by Hoivik et al. 
(2011) and Jowaed et al. (2010) in two separate studies 
where it was reported that intron methylation is associated 
with altered expression. Moreover, dense methylation sur- 
rounding transcription start site or near the first exon is 
tightly linked with gene silencing (Brenet et al. 2011). 

Till to date this is the first report of GMDS intron 
hypermethylation in chronic arsenic exposure with and 
without malignancy. Reports are also unavailable regarding 
association between p53, pl6 gene hypermethylation and 
GMDS gene hypermethylation in human cancer as well as 
in arsenic induced cancer. 

During the initial stage of the experiments we did 
observe some bands of hypomethylation, but we failed to 
clone them. It might be mentioned that in our previous 
investigations too, we observed far fewer hypomethylation 
cases. It is postulated that overexposure of arsenic and its 
biotransformation causes depletion of SAM, leading to 
hypomethylation of DNA. Hence extensive hypomethyla- 
tion probably needs a very high exposure, which is achieved 
in artificial tissue culture systems, but rarely in real life situ- 
ation. In the tissue culture experiments too, the study with 
cells exposed to arsenite for 2-4 weeks observed mostly 
hypermethylation and a few hypomethylation cases (Zhong 
et al. 2001). Chronic exposure of 18 weeks at low dose, on 
the other hand produced extensive hypomethylation and 
transformation in rat hepatocyte cell line (Zhao et al. 1997). 

Conclusion 

To sum up, this is the first report of GMDS gene 
fragment hypermethylation in the peripheral blood 
leukocyte DNA of persons exposed to arsenic. To ascertain 
this fragment of hypermethylation as a biomarker for arsenic 
induced cancer and chronic arsenic exposure researchers 
require repetition of such work in large sample group. 
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